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Abstract 

We derive a spectral representation for the two-point Green function for ar- 
bitrary composite field operators in Thermo Field Dynamics (TFD). A simple 



^ ■ way for calculating the spectral density within TFD is pointed out and com- 



pared with known results from the imaginary time formalism. The method is 
applied to hot 0^ theory. We give a compact derivation of the one- loop contri- 
bution to the shear viscosity and show that it is dominated by low-momentum 
plasmons. 
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Typeset using REVT[t;X 



I. INTRODUCTION 



In the context of heavy-ion colhsions at ultrarelativistic energies, transport properties of 
strongly interacting hot matter, in particular of the quark-gluon plasma (QGP), have become 
an issue of intense theoretical interest. The calculation of QGP transport coefficients from 
linear response functions through the Kubo-formalism ||T]J^ has so far been hampered by 
infrared problems associated with the non-linear selfinteractions among massless fields in 
hot QCD. A systematic and economic field theoretical approach, including resummation 
techniques for infrared divergences [^], to the calculation of linear response functions is 
clearly asked for. 

As an example for the type of objects one is interested in let us take the Kubo-formula 
for the shear viscosity 

V = -U d'x' f dt f dt' (n^'^iO), 7r,.(x', t')y^' = ^ lim^ 9p^j^,p) 

5 J J~oo J~oo 5 i^.P^O (juo 

Here n^j^^ is the traceless viscous-pressure tensor which is imbedded in the energy-momentum 
tensor and in the comoving frame has only spatial components, VTjj. For a scalar field 
T^ij = (SikSji — ^6ij6ki)dk(p di(j). Ptttt is the spectral momentum space density for the retarded 
thermal Green function of the composite field vr, as defined below. The efficient computation 
of this quantity will be the subject of this paper. 

The analytical structure of the retarded Green function {TT^'^{0),7[f^^{x' ,t')Y'^^ is similar 
to that of (0^(0), 0^(x', t'))^'^* p. The spectral density for the latter 2-point function was 
recently calculated at 1-loop order in the imaginary time formalism (ITF) by Jeon ^ . Here 
we repeat the calculation in the real-time formalism using the language of Thermo Field 
Dynamics (TFD) 0. Using the analyticity properties of the 2-point functions in ITF and 
TFD, we give a very compact derivation of the 1-loop contribution to the spectral density of 
the composite field 2-point function, Ptttt, in terms of the spectral density of the elementary 
field 2-point function, p^^. The latter was recently calculated at 2-loop order including 
infrared resummation Our approach follows closely the methods developed in PHTT 
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for correlation functions between elementary fields. It goes beyond the results reported in 
in that we include the full momentum dependence of the 2-loop scalar self energy. 



II. SPECTRAL REPRESENTATION OF COMPOSITE FIELD PROPAGATORS 

IN TFD 

For two arbitrary operators A, 13 the retarded thermal Green function is defined by 

D-j(ri -r2,ti = -^^(ti -t2)Tr(e^(^-^)[i(ri,ti),5(r2,t2)]±) , (2.1) 

where H is the hamiltonian, Q is the thermodynamic potential, and the plus (minus) sign 
applies for fermionic (bosonic) operators A, B. It can be expressed as an integral over the 
spectral density p^g{u!,p) in momentum space: 



DltiPo, P) = lini f ^du P^-^'^^P) . (2.2) 

ri^O+ J~oo Vn — UJ + iri 



Po - u + ir] 

In the Heisenberg picture the spectral density has the Lehmann representation 



p,,(c.,p) = (2vr)^Ee^^"~^"^ 

n,m 

{n\A{0)\m){m\B{0)\n) (1 - ae"^-) S{uj - E^^) 6'^'\p - p^„) , (2.3) 

where a = +1 (—1) for bosonic (fermionic) operators A, B, respectively. The sum is over 
a complete set of simultaneous eigenvectors of H and p, with eigenvalues En and p„, and 

Emn Eyyi Ejii Pmn Pm Prf 



The imaginary time Green function is obtained from Eq. (2J) by setting it = r. In 
momentum space it has a similar spectral representation, 



G,,(c.„,p)= r^du;^^^^, (2.4) 

J-oo lUn — UJ 



in terms of the same spectral function p^^. The periodicity of the ITF Green functions in 
imaginary time results in the discrete Matsubara frequencies Un = 2mxT and (2n+ l)7rT for 
bosonic and fermionic operators A and -B, respectively. The retarded Green function 
and the ITF Green function ( p.4|) are related by analytical continuation iuon ^ Po + if]'- 



^r«bo,p) = G',,K,p) 



iljJn^PO+iV 



(2.5) 



The propagators in TFD have a more comphcated structure: any physical field operator A 
is associated with a thermal ghost field operator A, with which it is combined into a doublet: 



'A' 
.Al 



(2.6) 



The TFD Green function is thus a 2 x 2 matrix 



Af^ (ri - r2, ti - h) = -I (0, /5|Ti[i"(ri, ti), B\v2, t2)]±|0, /3) , a, 6 = 1, 2 . (2.7) 



In momentum space it has the following spectral representation 
Afrpo,p)=/ du 



7,1 (^,P) ^7,'^(^>-p) 



(2.8) 



- + po + uj -tri_ 

where Tt is time ordering and the spectral function has the generalized Lehmann represen- 
tation 



7;^(^,p)=(27r)3^(0,/3|A'^(0)|n,m)(m,n|i?^(0)|0,/3)5(^-E„J5(3)^p_p^ 



, (2.9) 

7^^(a;,-p)^27r)3^(0,/3|5^(0)|n,m)(m,n|A'^(0)|0,/3)5(a;-E„J5(3)^p^p^j_ (2.10) 

n,m 

Here the thermal vacuum state |0,/3) is defined as 

|0,/5)=5]e'^(^-^^)/2|/t,^), (2.11) 
fe 

where fc) = |A;) ® and \k) is eigenstate of H with the same eigenvalue E^. 
Using Eq. (|2.11| ) and the 'tilde substitution law' 

{0,P\A\0)\n,m) = (-l)^^^^^-!)/^ exp[-/3E„„/2](0, /3|i(0)|n, m) , (2.12) 

{m,n\A\0)\0,P) = (-l)^^(^^+^)/2 gxp[-/3E„^/2](m, n|i(0)|0, /3) , (2.13) 

where is the fermion number associated with the operator A, and noticing that the 
physical field A only acts on the physical states \n) while the ghost field A acts only on the 
tilde states |n), 



Ck,l\A\n,m) = Skmm\n), (2.14) 
{k, l\A\n, fh) = 6in {k\A\Tfi) , (2.15) 

we can transform 7"^ in Eqs. ( p.9|) and ( p.lO| ) into similar form as Eq. ( |2.3| ). Some further 
algebra involving Eqs. ( p.2|) and ( |2.8|) then yields for the TFD propagator matrix the result 



\Apo,p)= (2.16) 
/ nret/ N 27ria exp(/3po/2) . . \ 

exp(/3po) - cr exp(/5po) - 

. s^sexp(/5po/2) . . / .N* 2m 

\ exp(/5po) - cr V AS / exp(/3po) - 0^ / 

where s^^ = (— l)^A{i^A+i)/2^ using Fa = —Fb since otherwise A^^ vanishes 0. 

With ( p^.2D this provides explicit expressions for each component of the 2-point Green 

function in terms of the spectral density p^g. The relationship of p^g to the (12)- and 

(21)-components of A^^ is particularly straightforward: 

Pab (Po, P) = ^ g/3po/2 ^AB (Po, P) ■ (2-17) 

Since the (12)-component of the 2-point Green function can be calculated directly using TFD 
Feynman rules, this relation provides a convenient way of obtaining the spectral density for 
the composite fields A, B. 

The above expressions apply in particular to the single-particle operators A = (f) and 
B = (j)^ {F^ = —Fg = 0, cr = 1). As shown in p^^^{uj, p) is real, depends only on |p| due 
to spatial isotropy, and is for any CPT-invariant equilibrium state an odd function of the 
frequency, 

(^,p) = e(^)p^^t(|^|,|p|), (2.18) 



where e{uj) = sgn(u;). Using Eq. ( p.5| ) and the following relation |T3[ for the analytically 
continued ITF propagator 

G^^^{po + iT],p) = G^^i{po + ie{po)r],p) + 2TTi9{-po)p^^t{\po\,p) , (2.19) 



we can write the full TFD propagator for the scalar field as 
A = A^^t(po,p) = 



(2.20) 



exp /? po - 1 



-2m 



-2Txi 



exp(/5bo|/2) 



exp(/?|po|/2) 
exp(/?bo|) - 1 "^"^ 

2vrz P00t(bo|,p) 



V "••^ exp(/jbo|)-l '^'^^'^'^°''^^' -G';0t(Po + ^e(PoKp)-- 
with the following representations in terms of the ITF self energy S: 
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exp(/3bo|)-l / 



jp. + ReS(bo|,p) +^ImS(bo|)P) 
.(bo|,p) = -^ 



G^^^ ( bo I + ^'7> P) - G^^^ ( bo I - P) 



ImS(bo 


,P) 


b^ — + Re S(bo 


,p)]' + [ImE(bo 


,P)f 



For a free scalar field both real and imaginary parts of S vanish, and the spectral 
,t degenerates to a (5-function: 



— lim 



Eq.( ^.2(J| ) in this case reduces to the free TFD propagator 

/ 1 2m5{p^-m'^) _27ri ^(p^ - m^) exp(/?bo|/2) \ 



(2.21) 

(2.22) 

density 

(2.23) 



Ao 



p2 _ ^2 _|_ exp(/3bo|) — 1 
2m5{p^ - m^) exp(/5bo|/2) 



exp(/?bo|) - 1 

27ri 6{p'^ — m^] 



V 



p2 _ ^2 



exp(/5bo|) - 1 

The full TFD propagator matrix A satisfies the Dyson equation [1^ 
(A-i)"'' = (p2-^V' + S'^'(Po,p), 



(2.24) 



ir] exp(/3bo|) - 1 / 



(2.25) 



where r = diag(l, —1). This defines a 2 x 2 TFD self energy matrix S . Comparing this 



with the inverse of the matrix ( |2.2CI|) , it is easy to obtain the following relations: 

=/3po _ 1 g/3po _ 1 



Im E 



Im S 



+ 1 ''^^ 2e/^Po/2 



ImS 



12 



Re S = Re S 
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-ReS 



22 



Re T}"^ = Re S^^ = . 



(2.26) 
(2.27) 
(2.28) 



They were previously obtained in [^|T5| by different methods. The authors of Ref. [|T5 
pointed out that the calculation of ImE from ImS^^, using TFD Feynman rules for the 
evaluation of the latter, is particularly convenient and showed explicitly that the result is 
identical to the evaluation in ITF. 

It is straightforward to generalize the expressions ( |2.2(j| ) and ( f^.24D as well as the self 
energy relationships (|2.26| - |2.28| ) to the case of fermionic single particle operators. 



III. SHEAR VISCOSITY IN HOT ch-^ THEORY 



We will now use these results to give a five line derivation of the spectral function for the 
composite field 0^ within hot theory. We concentrate on the same quantity as studied 
in Ref. 0, namely 



T]^2^2 = lim 

p,po^O 



Po 



(3.1) 



whose relation to the shear viscosity (|1 . 1|) was mentioned above. From Eq. ( p.l7| ) we obtain 

i e^Po - 1 



Ai2.2(po,P) 



(3.2) 
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The lowest order contribution to Ap^2(po)P) is shown diagrammatically in Fig. 1. This 
is a skeleton diagram, i.e. the full single-particle propagators ( |2.2(J| ) must be used for the 
internal lines Jl],^. For the one-loop diagram in Fig. 1 only its (12)-component is needed. 
We find 



i A^2^2(po,p) = 2 



exp(^)p^^(A;o,k) 
exp{pko) - 1 



exp(fcM; 



P^<t>iPo -ko,p- k) 



exp(/?(po - ko)) - 1 



2e~ 



with the Bose distribution 



fih) f{po - h) P<t><t>{ko, k) p4,<i,{po - fco, P - k) , (3.3) 



f{ko 



1 



Combining this with Eqs. 



exp(/3fco) - 1 
and (B.ll) we obtain 



(3.4) 



where we used Eq. (|2.18| ) as well as the identity /(— /co) = ~[1 + /(^o)]- The factor /? results 
from lim (e^^" — l)/po- Up to a factor of 27r resulting from our different normalization ( p.22| ) 

po^O 

of the spectral density, this result is identical with the (much lengthier) calculation ||^ from 
the cutting rules in ITF. 

To further evaluate Eq. ( p.5| ) we need the spectral density p^(f,{kQ, k) for the full single 
particle propagator. Since the calculation of ^1^2^°^ requires taking the zero momentum 
limit of the loop diagram, both internal lines in Fig. 1 can become soft. For massless A^"^ 
theory this means that, in order to avoid infrared divergences, we should use resummed 
effective propagators 0,0. The so-called resummation of "hard thermal loops" |Q in this 
case generates a thermal mass for the scalar field which acts as an infrared cutoff. Using 
such a resummation scheme, we recently performed |§] a 2-loop calculation of p^^. We found 
that for weak coupling among the scalar fields (A/24 ^ 1) the spectral function is sharply 
peaked around the plasmon frequency 



a;,(k) = y'k^ + - ^ ^ V^k^T^ (3.6) 

where 

mth = (3.7) 

is the "thermal mass" and rrip is the plasmon mass 0. The spectral function p,^^ can then 
be expressed to good approximation in the form of a relativistic Breit-Wigner function, 

1 2fco7(k) 

p^^ih, k) ^ - (^2_^2(k))2+4^2^2(k) ' (3.8) 

where 

= i'-^M^m (3.9) 

2u;p(k) 

is the on-shell damping rate for the scalar plasmon. Its leading contribution comes from the 
2-loop diagram in Fig. 2 and is given by 
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7(k) 



1 /-Ikl 



256 7r3a;p(k) |k| Jo 



dq 



' i-C 



u 



' (e-c)(i-ec) ' 



where 



and Ij2{z) is the Spence function 

u{z) ^ - r 

io 



ln(l - t) 



Substituting Eq. ( p.8|) into Eq. ( p.5| ) and integrating over fco we obtain [p],|5|,p!7| 

l-loop /K(k)) [1 + /K(k))] 

^-^'^^ ^^2717 (27r)3 cu2(k)7(k) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



IV. RESULTS AND DISCUSSION 



We will now evaluate the momentum integral in ( |3.13| ). It can be studied analytically 
in the weak coupling limit A ^ 1. Let us first consider an approximation which has been 
investigated before where the momentum dependence of the plasmon width is neglected 
and 7(k) is replaced by its zero-momentum limit 



7(0) 



Im S(mp, 0) 



A2 r2 



2mr 



256 7r3 



A^T^ 



Inserting this into Eq. ( |3.13|) and integrating by parts leads to 

384 dx 1 



1536'7rmr 



l + 0(v^lnA)l . (4.1 



1— loop 



(4.2) 



7{k)=7{0) T A2 Jo (x2 + a^)^/^ e^^^ - 1 ' 
with a = rrip/T. For small a <^ 1 the integral can be evaluated following the methods of 



Appendix C in Ref. |T8[, and we obtain W9 



1— loop 



_ 96 1 

r{k)=7{0) ^ 



677" 



+ 0{X)\ . 



(4.3) 



With the exphcit momentum dependence |^ of the plasmon decay width from Eq. ( p.lO| ) 
this approximation can be avoided. Inserting Eq. ( |3.10[ ) into Eq. p.l3| ) and integrating by 
parts, one finds 

i-ioop ^ dx A{x;a) - xA'{x;a) 

where x = k/T, the prime denotes d/dx, and 

256 TT^ 128 TT^ 

A{x;a) = ^^ujp{k)^{k) = Im S(o^p(k), k) (4.5) 



- / dz 



X Jo 

with (see Eq. ( PID ) 



^_g-v^^^ ^_^_v^^^ a = ^. (4.6) 

The function A{x\ a) is shown in Fig. 3. For large momenta, x oo, it approaches the 
constant value |16| 

— 0,^ de 



lim Aix- a) = - r dz In (l- e'^^^) = f 

Jo ^ ^ Ja 

7r^/6 for a — * ; 
aKi (a) for a ^ 1 . 



- 1 



(4.7) 



At zero momentum it takes the value 16 



A(0;a) =L2(e-'^) (4. 



which again approaches in the weak coupling limit a=yA/24^0. Unfortunately, in 
the weak coupling limit the function A{x; a) exhibits a very strong momentum dependence 
in the region a < x < 1; for a = the limit x — is non-analytic. If we nonetheless neglect 
the momentum dependence of A, Eq. ( [4.4| ) becomes very simple: 

i_ioop _ _6£ _6_ /•°° dx 

In the limit a ^ this integral diverges logarithmically in the infrared. The singular 
behaviour can be isolated by writing 
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dx /■°° e de r°° de e - \/e^ - de 



ev"^^ - 1 J a ^Je^ - a2 - 1 7a - 1 7a - - 1 

= lnf-')+ln2 + 0(a) + of^) . (4.10) 



In the second line we cut off the second integral at the upper end at a point 5 with a <^ 
5 <^1. This allows to approximate the Bose distribution as 1/e. One easily verifies that the 
remaining integral from 5 to oo is finite and of order as indicated. We thus find with 

this approximation [|T^ 

384 1 A / T \ . , „ „^ /m, ■ 



192 1 



^ fin (\] + ln(96) + O i^fx)] . (4.11 



7r2rA2V VA. 

The additional logarithmic divergence in the weak coupling limit of ( [4.11 ) compared to 



is generic; its coefficient and the next-to-leading constant term depend, however, on 
the momentum dependence of A{x\a). Due to the non-analytic behaviour of A(x; 0) near 
a; = we were not successful in extracting an analytical expression similar to ( [1.1 1| ) for the 
full integral (|4.4| ). The numerical result shown by the dots in Figure 4 can, however, be 
excellently fit by 



1— loop / 1— loop 



0.3282 In ( — ) + 1.41682 ^ - In ( — ) + ^2 (4.12) 



7{k)=7{o) \ nir, } TT \m 

which suggests the analytical behaviour 



1— loop 



96 1 / / T \ /- _ /m 



,2^2 



Tc'^T A2 y \^p / V T 

4J 



^ fin (\] + ln(24) + 27rv^ + O (\^)] . (4.13) 



7r2TA2V VA. 

As intuitively expected, the viscosity decreases as the coupling strength increases: In 
relaxation time approximation all transport coefficients are proportional to the relaxation 
time, which again is inversely proportional to the scattering rate which grows with the 
coupling strength. On the other hand, the rate at which our 77,^2^2 decreases with the 
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coupling strength A is different from the behaviour of the physical shear viscosity rj. As 
discussed in the Introduction, the latter involves four additional spatial derivatives acting 
on the scalar field. This translates |1| into an additional factor k'^ in the integrand of 
Eq. ( |3.13| ), which removes the infrared divergence of this integral in the limit a — > and 
thereby also the leading logarithmic term in our final result ( [4.13| ). This observation provides 
a partial explanation for the qualitatively different behaviour of the viscosity as a function 
of mp/T which was recently found by Jeon in a more complete study of g(j)^ + theory 
iBf : his numerical results indicate an increase of X^rj with increasing A. 

Still, as first pointed out by Jeon in [0] and then quantitatively analyzed in |16| , the fact 



that the shear viscosity rj is proportional to 1/A^ raises a serious problem: simple power 
counting arguments show that then all planar ladder diagrams of the type shown in 
Fig. 5 can also contribute to the leading order result for the viscosity. The summation 
of this infinite series of ladder diagrams is nontrivial and has been recently performed by 
Jeon [|1^] using the imaginary time formalism. This work was a genuine tour de force, and 
the prospect of generalizing it to non-abelian gauge theories like QCD seems frightening. 
However, based on the results presented in this work, we believe that the methods developed 
here in the framework of TFD will help to redo the analysis in a more efficient way and allow 
for an easier generalization to QCD. Work in this direction is in progress. 
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FIGURES 



(ko,k) 




iU-ko'P-k) 

FIG. 1. One-loop skeleton diagram for A^^2{po,p)- The heavy lines denote full single particle 
propagators. 




FIG. 2. The leading contribution to the imaginary part of self energy ImS. The dotted lines 
denote resummed effective propagators. 
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FIG. 3. (a) The function {<o/it'^)A{x\ a) vs. a; for a = 0.01, 0.1, and 0.4 (from top to bottom), 
(b) Close-up of (a) for small x. 




FIG. 4. The dots represent numerical results for the ratio between V^2^2°^ and 77^2^°°^ 
as a function of rUp/T. The solid line indicates a fit of the form (I/tt) ln(T/mp) + \/2. 



7(k)=7(0) 
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